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ABSTRACT 

Aims. The rapidly rotating Be star p Persei was spun up by mass and angular momentum transfer from a now stripped-down, hot 
subdwarf companion. Here we present the first high angular resolution images of p Persei made possible by new capabilities in long- 
baseline interferometry at near-IR and visible wavelengths. We analyzed these images to search for the companion, to determine the 
binary orbit, stellar masses, and fluxes, and to examine the geometrical and kinematical properties of the outflowing disk surrounding 
the Be star. 

Methods. We observed p Persei with the MIRC and VEGA instruments of the CHARA Array. MIRC was operated in six-telescope 
mode, whereas VEGA was configured in four-telescope mode with a change of quadruplets of telescopes during two nights to improve 
the (u, v) plane coverage. Additional MIRC-only observations were performed to track the orbital motion of the companion, and these 
were fit together with new and existing radial velocity measurements of both stars to derive the complete orbital elements and distance. 
We also used the MIRC data to reconstruct an image of the Be disk in the near-IR H- band. VEGA visible broadband and spectro- 
interferometric Her observations were fit with analytical models for the Be star and disk, and image reconstruction was performed on 
the spectrally resolved Her emission line data. 

Results. The hot subdwarf companion is clearly detected in the near-IR data at each epoch of observation with a flux contribution of 
1.5% in the H band, and restricted fits indicate that its flux contribution rises to 3.3% in the visible. A new binary orbital solution is 
determined by combining the astrometric and radial velocity measurements. The derived stellar masses are 9.6±0.3M o and 1.2±0.2M o 
for the Be primary and subdwarf secondary, respectively. The inferred distance (186 + 3 pc), kinematical properties, and evolutionary 
state are consistent with membership of p Persei in the a Per cluster. From the cluster age we deduce significant constraints on the 
initial masses and evolutionary mass transfer processes that transformed the p Persei binary system. The interferometric data place 
strong constraints on the Be disk elongation, orientation, and kinematics, and the disk angular momentum vector is coaligned with 
and has the same sense of rotation as the orbital angular momentum vector. The VEGA visible continuum data indicate an elongated 
shape for the Be star itself, due to the combined effects of rapid rotation, partial obscuration of the photosphere by the circumstellar 
disk, and flux from the bright inner disk. 

Key words, stars: individual: p Persei - stars: binaries 


1. Introduction 

For about 30 years, Be stars have been favorite targets of op¬ 
tical interferometers (Thom et al. 1986; Mourard et al. 1989), 
and they are important objects to study the role of rotation in 
the link between a star and its close environment. Be stars are 
rapidly rotating, main-sequence B-stars that eject gas into cir¬ 
cumstellar disks (Rivinius et al. 2013). The “e” suffix in their 

* Based on observations with MIRC-6T and VEGA-4T instruments 
on the CHARA Array 


spectral classification refers to the presence of Balmer emission 
lines that form in their outflowing disks. Such emission appears 
in B-stars with equatorial velocities above approximately 75% 
of the critical rate, and processes related to nonradial pulsation 
and/or small-scale magnetic fields probably aid mass loss into 
the disk (Rivinius et al. 2013). In some circumstances, these 
processes may produce high-energy X-rays. The Be star y Cas- 
siopeiae, for example, exhibits thermal X-ray emission. In a re¬ 
cent series of two papers on y Cas, Smith et al. (2012) and Stee 
et al. (2012) have reviewed the most remarkable features of this 
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star: X-ray activity, magnetic field, critical rotation, Keplerian 
rotation, disk elongation, and binarity. These papers have also 
shown once more that a multiwavelength and multitechnique ap¬ 
proach is key to making progress in understanding the physics of 
these complex systems. 

The origin of the rapid spin of Be stars is unknown. Most Be 
stars appear to be somewhat evolved objects (McSwain & Gies 
2005; Zorec et al. 2005), so that their fast spin represents pro¬ 
cesses that occur well into their main-sequence lives. One pos¬ 
sibility is that B-stars spin up as they conclude core H-burning 
due to the redistribution of internal angular momentum (Ekstrom 
et al. 2008; Granada et al. 2013; Granada & Haemmerle 2014). 
A second possibility is that Be stars were spun up by mass and 
angular momentum transfer in a binary (Pols et al. 1991; de 
Mink et al. 2013). In this case, the companion would lose most 
of its envelope and would appear as a hot, stripped-down He-star 
remnant. The Be binary star ip Persei represents the first detec¬ 
tion of a Be star with a hot, faint companion (Poeckert 1981; 
Thaller et al. 1995; Gies et al. 1998), and it is one of three such 
systems known at present (Peters et al. 2013). Clearly it is im¬ 
portant to study such Be binaries to determine the physical prop¬ 
erties of the stars at the conclusion of their transformative inter¬ 
action. 

This paper is an analysis of the results obtained on ip Persei 
thanks to a combined visible and infrared interferometric cam¬ 
paign with the instruments VEGA (Mourard et al. 2009) and 
MIRC (Monnier et al. 2008) on the CHARA Array (ten Brum- 
melaar et al. 2005) and is complemented with a new analysis of 
a collection of radial velocity measurements. In Sect. 2 we sum¬ 
marize the known facts on ip Persei relevant to this study, and we 
compile radial velocity measurements in Sect. 3. The journal of 
interferometric observations is presented in Sect. 4, and MIRC 
data are presented in Sect. 5 along with a global astrometric and 
radial velocity analysis. The VEGA data and analysis are pre¬ 
sented in Sect. 6. The scientific implications of our results are 
presented and discussed in Sect. 7, and we conclude in Sect. 8. 

2. Main characteristics of ip Persei 

The source ip Persei (HD10516, HIP8068) is a bright star (m v = 
4.09) in the northern sky that is classified as B 1.5 V:e-shell (Slet- 
tebak 1982). The Hipparcos parallax (van Leeuwen 2007) is 
7i - 4.54 ± 0.2 millisecond of arc (mas), which gives a distance 
d = 220 ± 10 pc. Using the interstellar Ca II distance scale, 
Megier et al. (2009) estimate a distance of d = 195 + 32 pc. 

Campbell (1895) published a list of 32 stars presenting both 
bright and dark hydrogen lines and addressed the question of 
explaining these interesting and varying features, ip Persei ap¬ 
pears in this list, and about 120 years later, we are still consid¬ 
ering it as a puzzle. Until 1960, most of the papers published 
on ip Persei after this discovery were dedicated to samples of 
stars and statistical studies. More detailed work followed this 
period, and Slettebak (1966) (confirmed later on by Abt et al. 
2002) published a measurement of the projected rotational ve¬ 
locity of V sin/ = 410-450 km s -1 , which is estimated to be 
around 85% of the equatorial breakup velocity. 

Poeckert (1981) presented the first evidence of a compan¬ 
ion in ip Persei and identified it as a helium star. Based on 
the velocity curves, he estimated masses of M p = 21 M© and 
M s = 3.4 M q . The primary is the Be star, whereas the sec¬ 
ondary is the hot companion with an effective temperature es¬ 
timated to be around 50000K in order to ionize helium. The 
photospheric spectrum of the hot companion was first detected 
in spectra from the International Ultraviolet Explorer satellite 


by Thaller et al. (1995), who showed that the spectrum appears 
similar to that of the sd06 subdwarf star HD 49798. Bozic 
et al. (1995) presented a comprehensive review of the avail¬ 
able photometric and radial velocity data, and they proposed an 
orbital solution with a period P = 126.6731 days, an epoch 
of Be star superior conjunction T sc = HJD2435046.73 (HJD 
= heliocentric Julian Date), and revised masses: M p between 
16 and 22 M© and M s between 1.7 and 2.2 M©. They clas¬ 
sified the primary star as B0.5e. The companion was finally 
confirmed by Gies et al. (1998) thanks to high-resolution UV 
spectroscopy with the Hubble Space Telescope. The subdwarf- 
to-Be flux ratio is 0.165 ± 0.006 (resp. 0.154 ± 0.009) for the 
1374A region (resp. 1647A). Moreover, they derived a double- 
lined solution for the radial velocity curve that yields masses of 
M p = 9.3 ± 0.3 Mq and M s = 1.14 ± 0.04 M©. The lower mass 
estimates with respect to the results from Bozic et al. (1995) are 
explained by a much smaller secondary semiamplitude in the 
RV measurements, which leads to a smaller projected semima¬ 
jor axis and thus lower masses. The subdwarf effective tem¬ 
perature is T q ff = 53000 ± 3000X and the surface gravity is 
logg = 4.2 ± 0.1. 

Through the study of the helium lines, Stefl et al. (2000) ar¬ 
gued for an origin of the emission in the outer parts of the disk 
surrounding the primary star. Their observations agree with a 
scenario where the outer parts of an axisymmetric disk are illu¬ 
minated by the radiation of the secondary. They also favor in¬ 
homogeneities in the global density pattern of the inner regions 
as evidenced by the emission line asymmetry. The model of the 
FeII5317A and He I 6678A and 5876A emission lines of ip 
Persei proposed by Hummel & Stefl (2001) yields the size and 
shape of the excitation region in the circumprimary disk. The 
Fe II emission originates within 12 stellar radii in an axisym¬ 
metric disk around the primary, whereas the He I emission is best 
fit by a disk sector with a radius of 15 stellar radii and opening 
angle of - 120° facing the secondary (Hummel & Stefl 2003). 
They also confirmed that the complex structure of the emission 
line profiles is due to an external illumination. 

Polarimetric studies of ip Persei were made in the late 1990s 
through the work of Clarke & Bjorkman (1998) and Ghosh et al. 
(1999). A wavelength dependence of the polarimetric fluctua¬ 
tions was detected and associated with events deep inside the 
binary system, the amount of scattered radiation being primarily 
controlled by the opacity of the primary star’s disk. Ghosh et al. 
(1999) measured an intrinsic polarization of ip Persei equal to 
1.57+0.11% (resp. 1.36+0.10%, 1.19+0.09% and 0.98+0.08%) 
in the B (resp., V, R and I) band. The position angle (measured 
east from the north celestial pole) of the scattering disk normal 
axis is PA = 27 + 2°. 

Optical interferometric observations of ip Persei (see Ta¬ 
ble 1) were first reported in 1997 with the Marklll interferom¬ 
eter (Quirrenbach et al. 1997). Using a wide spectral band of 
54A centered on H a, they fitted their data with a model of a 
Gaussian disk plus a stellar uniform disk of 0.39 mas contribut¬ 
ing 23% of the flux in the spectral band. Using the NPOI in¬ 
terferometer with an H a filter of width 39A+3A , Tycner et al. 
(2006) used the same model as Quirrenbach et al. (1997) and 
obtained similar results. The authors also derived an inclination 
angle i > 74°. Finally, Gies et al. (2007), using the CHARA 
Array in the X'-band, fitted a model with a Gaussian elliptical 
disk plus a central source and a secondary object accounting for 
6% of the flux (A K' = 2.9 mag). The flux contribution from the 
central source is estimated as 53.8 + 1.5%. 
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Table 1 . Interferometric estimates of disk parameters. References are 
(1) Quirrenbach et al. (1997), (2) Tycner et al. (2006), (3) Gies et al. 
(2007), (4) Touhami et al. (2013), (5) this paper. Columns 1 and 2 list 
the ratio of the minor to major axis r and its uncertainty, Columns 3 and 
4 list the disk position angle along the major axis PA and its uncertainty, 
and Cols. 5 and 6 list the angular FWHM of the disk major axis 0 m j and 
its uncertainty. 


r 

Sr 

PA 

SPA 

$maj 

<S0maj 

Band 

Ref. 

0.46 

0.04 

-62 

5 

2.67 

0.20 

Ha 

a) 

0.27 

0.01 

-61 

1 

2.89 

0.09 

Ha 

(2) 

0.00 

0.22 

-44 

3 

2.30 

0.08 

K' 

(3) 

0.10 


-44 

4 

2.44 

0.22 

K f 

(4) 

0.15 

0.09 

-64 

3 

1.44 

0.41 

H 

(5) 

0.37 

0.03 

-76 

1 

0.48 

0.05 

R 

(5) 

0.34 

0.03 

-69 

2 

0.82 

0.10 

Ha 

(5) 


The main features of the system may be summarized as fol¬ 
lows: The system is a known spectroscopic binary with a period 
of almost 127 days. The primary object is a Be star with a very 
high rotational velocity (450 km s -1 ) that presumably presents 
a very elongated photosphere seen almost equator-on. As ev¬ 
idenced by the first interferometric measurement and in accor¬ 
dance to its shell classification, it is expected that the inclination 
angle is at almost 90°. The companion is a hot subdwarf star. 
We can estimate from the UV flux a contribution of 8% or less 
in the V - and H- band continuum. The main goals of our study 
are to attempt a direct detection of the companion in the optical 
to determine its physical and orbital properties and to continue 
the study of the disk geometry and kinematics. 

3. Spectroscopic orbital determination from radial 
velocity measurements 

A combined analysis of the astrometric and radial velocity orbits 
provides the key step to determining the stellar masses and dis¬ 
tance. However, a comprehensive spectroscopic orbit was last 
determined by Bozic et al. (1995), so it is worthwhile to re¬ 
assess the spectroscopic orbital elements based on observations 
that are more contemporaneous with the interferometric obser¬ 
vations. Here we present new radial velocity measurements for 
the Be primary and for the elusive sdO secondary. 

Bozic et al. (1995) showed that the wings of the H a emission 
line that are formed in the inner disk close to the Be star display 
radial velocity variations of the Be star. Thus, we collected the 
available H a spectroscopy from the Be Star Spectra Database 1 
(Neiner et al. 2011) and from the work of Grundstrom (2007). 
These were supplemented by several spectra we obtained with 
the KPNO Coude Feed telescope and one spectrum published 
by Tomasella et al. (2010). All of these are high S/N spectra 
with a resolving power generally better than R = A/AA = 5000. 
The reduced heliocentric Julian dates (RJD) of observation and 
origins of these 79 spectra are listed in Table 2. The H a emis¬ 
sion strength changed by a factor of two during these observa¬ 
tions (with local maximum strength in 2005 and local minimum 
strength in 2011), and there were changes visible in the shape 
of the line core. However, the shapes of the emission wings re¬ 
mained more or less constant throughout this period, and this 
suggests that fast-moving gas close to the base of the disk at the 
photosphere was present throughout. Consequently, measure¬ 
ments of the Doppler shifts from the wing bisector should pro¬ 
vide a reasonable estimate of the motion of the gas closest to the 

1 http://basebe.obspm.fr/basebe/ 


Be star itself. Bozic et al. (1995) came to a similar conclusion 
and found no evidence of systematic deviations in the resulting 
velocities (which might be present, for example, if high-speed 
gas streams contributed to the emission flux). 

We measured the H a wing velocity by determining the emis¬ 
sion line bisector position using the method described by Shafter 
et al. (1986). The radial velocity is found by cross-correlating 
the observed profile with a pair of oppositely signed Gaussian 
functions offset to positions of ±250 km s _1 from line center and 
then measuring the zero-crossing position of the resulting cross¬ 
correlation function. The uncertainties in the resulting measure¬ 
ments were estimated using the semi-empirical expression from 
Grundstrom (2007), but we caution that these estimates do not 
account for long-term changes in the emission profile that also 
influence the wing measurements. In most cases we were able 
to also measure accurate positions of the nearby atmospheric tel¬ 
luric lines, and we used these to make small corrections to the 
wavelength calibration in the topocentric frame. The resulting 
radial velocity measurements appear in Col. 3 of Table 2. For 
convenience, we also include a set of six radial velocities mea¬ 
sured by Gies et al. (1998) for the Be star photospheric lines in 
far-ultraviolet spectra obtained with the Hubble Space Telescope. 

We made a single-lined spectroscopic orbital determination 
from these radial velocities alone using the program SBCM 
described by Morbey & Brosterhus (1974). We began mak¬ 
ing circular fits with the period derived by Bozic et al. (1995), 
P = 126.6731 ± 0.0071 d, and we derived an orbital semiampli¬ 
tude of K a = 11.4 ±1.0 km s -1 . This agrees within the uncertain¬ 
ties with the value reported by Gies et al. (1998), K a = 10.0 + 0.8 
km s" 1 . 

The measurement of radial velocities for the sdO companion 
is much more difficult. There are several features in the optical 
spectrum that display the antiphase Doppler shifts expected for 
the companion, but these are usually associated with gas near 
the companion star. These include weak He II 4686A emis¬ 
sion that may form in a disk close to the companion (Poeckert 
1981), He I 4026A , 447lA “shell” lines formed farther away 
(Poeckert 1981), and the He I 6678A emission line that proba¬ 
bly originates in the outer disk of the Be star facing the hot sub¬ 
dwarf (Stefl et al. 2000; Hummel & Stefl 2001). There are only 
six HST measurements of the photospheric lines of the compan¬ 
ion itself that were made in the far-ultraviolet spectrum by Gies 
et al. (1998), and we include these in Table 2. There are also 16 
high-resolution FUV spectra of cp Persei available in the archive 
of the International Ultraviolet Explorer Satellite (Thaller et al. 
1995), and we successfully detected the spectral signature of the 
companion in each of these using the cross-correlation methods 
described by Peters et al. (2013). We measured radial veloci¬ 
ties for the companion using the same model template spectrum 
adopted in the work of Gies et al. (1998), and our derived radial 
velocities appear at the bottom of Table 2. The uncertainties as¬ 
sociated with these velocities are probably larger than reported in 
the table. Indeed, since they are too faint, the interstellar lines in 
the spectrum of cp Persei cannot be used to correct for zero-point 
shifts of the star position with respect to the slit in the large- 
aperture observations (Peters et al. 2013). 

We made a single-lined spectroscopic orbital solution for the 
secondary again using the orbital period from Bozic et al. (1995), 
and we obtained a semiamplitude of K b = 83.6 ±2.1 km s _1 that 
agrees within errors with the determination from the HST spec¬ 
tra of K b = 81.2 ± 0.6 km s _1 (Gies et al. 1998). Elliptical fits 
provided no improvement in the residuals. Thus, the IUE mea¬ 
surements tend to confirm the results from the HST spectra. On 
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the other hand, because the measurement errors are larger for the 
IUE velocities, it is not clear whether or not including them im¬ 
proves our estimate of the semiamplitude K h . In the following 
section, we restrict our considerations to the HST velocities of 
the companion. Note that the FUV derived estimate of the sec¬ 
ondary semiamplitude is somewhat lower than that found from 
optical studies (K h = 105 ± 7 km s -1 from He II 4686A , Poeck- 
ert 1981; Kb « 90 km s -1 from He I 4026A , 4471A , Poeckert 
1981; K b = 101 ± 6 km s" 1 from He I 6678A , Bozic et al. 1995) 
presumably because the latter are influenced by gas motions and 
a location of origin that differ from that of the hot companion. 
The He I 6678A emission, for example, probably originates in 
the outer disk facing the companion at a position closer to the Be 
star and thus with a higher Keplerian orbital speed than that of 
the companion. 

4. Interferometric observations and data analysis 
principles 

Interferometric observations were carried out on the CHARA 
Array (ten Brummelaar et al. 2005), located at Mount Wilson 
Observatory. It consists of six 1 m telescopes sending their 
beams, after compression, through vacuum pipes to the recom¬ 
bining lab, through the fixed and movable delay lines with the 
correction of longitudinal dispersion, and the beams are re¬ 
arranged before feeding the instruments. For the purpose of this 
work we have used the MIRC instrument (Monnier et al. 2004, 
2008) operated in the H band using all six CHARA telescopes 
(Che et al. 2012a) and the VEGA spectrograph used in the 4-T 
mode (Mourard et al. 2011). The spectral configurations used 
with the different instruments are summarized in Table 3. The 
different configurations we used during these observations are 
presented in Table 4. The corresponding ( u , v) coverage, gener¬ 
ated with the Aspro2 service 2 , is presented in Fig. 1. 

Table 3. Spectral configurations used with VEGA and MIRC for the 
<p Persei observations. Column 2 gives the spectral band, Col. 3 the 
spectral width of one channel, and Col. 4 the number of channels. 



A nm 

AT nm 

Channels 

VEGA 

[483;498] 

15 

1 


[636;651] 

15 

1 


[648;663] 

15 

1 


[663;678] 

15 

1 

VEGA 

[646;666] 

0.1 

200 

MIRC 

[1501; 1712] 

35 

6 


For the VEGA observations, VEGA and MIRC were used 
simultaneously, MIRC with six telescopes in H band and VEGA 
with two different four-telescope configurations. MIRC was 
used as a 6-T group delay sensor, and corrections were sent di¬ 
rectly to the main delay lines. An initial cophasing of the two 
instruments was made using a reference target thanks to the fine 
adjustments of the MIRC feeding optics and corresponding cor¬ 
rections on the main delay lines. The observing sequence was 
then 1) MIRC finds and stabilizes the fringes, 2) MIRC and 
VEGA start recording data, 3) at the end of the VEGA sequence, 
MIRC does its shutter and photometry sequence. 

Two different calibrators were used for the VEGA obser¬ 
vations: HD3360 and HD25642. They were identified us¬ 
ing the SearchCal service 3 (Bonneau et al. 2006, 2011) and 

2 Available at http://www.jmmc.fr/aspro 

3 Available at http://www.jmmc.fr/searchcal 
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Fig. 1. Top: (w,v) plane coverage for MIRC. Bottom: (w,v) plane 
coverage for VEGA. The two figures use the same spatial frequency 
scale. 

the CDS Astronomical Databases SIMBAD and VIZIER 4 . We 
used the SearchCal spectrophotometric uniform disk angular 
diameter determinations in the R and H bands, which led to 
Qjjd-r — 0.284 ± 0.020 mas and 0jjd-h ~ 0.287 ± 0.020 mas 
for HD3360 and to 0ud-r = 0.439 ± 0.031 mas and 0 UD - H = 
0.447 ± 0.031 mas for HD25642. In addition to these calibra¬ 
tors, MIRC calibration employed the additional calibrators 7 
And (0.65 ± 0.03 mas), 6 Cas (0.57 ± 0.04 mas), and HD 33167 
(0.50 ± 0.04 mas). 

The data processing followed the standard MIRC (Monnier 
et al. 2007, 2012) and VEGA (Mourard et al. 2009, 2011) pro¬ 
cedures. Based upon our understanding of the current VEGA 
limitations (Mourard et al. 2012), we decided to reject all clo¬ 
sure phase measurements with VEGA and to set a conservative 
lower limit to the absolute uncertainty of the squared visibility 
measurements at the level of 0.05. 

5. Near-IR image reconstruction and companion 
detection 

5.1. MIRC results 

For the purposes of this paper, we wished to extract two major 
findings from the MIRC data. First, we sought to directly de¬ 
tect the companion for the first time to determine the visual orbit 
of the (f Persei binary system. Secondly, we aimed to take ad- 

4 Available at http://cdsweb.u-strasbg.fr/ 
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Fig. 2. (left) Here we show the two-dimension visibility function for ip Persei from MIRC. We combined all our epochs together and averaged 
with a 30m smoothing window, maintaining each wavelength channel, (right) This shows the same averaged data as a function of radial spatial 
frequency unit. In addition, we have included model visibilities for position angles along the major and minor axes of the average disk model (see 
Table 5). 


Table 4. ip Persei observations with VEGA and MIRC on the CHARA 
Array. Column 1 gives the UT date, Column 2 the reduced heliocentric 
Julian date RJD, Col. 3 the orbital phase of the companion according 
to our newly determined ephemeris, Col. 4 the telescope configuration, 
and Col. 5 the UT time of observation. 


Date 

RJD 

Phase 

Telescopes 

UT 

VEGA 

2011Sep28 

55832.86 

0.81 

S2E2W1W2 

08:40 

2011 Oct 18 

55852.77 

0.97 

S1S2W1W2 

06:26 

2011 Oct 18 

55852.85 

0.97 

S1S2W1W2 

08:31 

2011 Oct 18 

55852.96 

0.97 

E1E2W1W2 

10:58 

2011Octl9 

55853.72 

0.98 

S1S2W1W2 

05:16 

2011Octl9 

55853.78 

0.98 

S1S2W1W2 

06:39 

2011 Oct 19 

55853.88 

0.98 

E1E2W1W2 

09:02 

2011Octl9 

55853.93 

0.98 

E1E2W1W2 

10:21 

MIRC 

2011Sep03 

55807.91 

0.62 

S1S2E1E2W1W2 

09:50 

2011Sep28 

55832.86 

0.81 

S1S2E1E2W1W2 

08:38 

2011 Oct 18 

55852.87 

0.97 

S1S2E1E2W1W2 

08:52 

2011Octl9 

55853.76 

0.98 

S1S2E1E2W1W2 

06:14 

2012Augl8 

56158.02 

0.38 

S1S2E1E2W1W2 

12:29 

2012Sepl5 

56185.81 

0.60 

S1S2E1E2W1W2 

07:26 

2012Oct31 

56231.77 

0.96 

S1S2E1E2W1W2 

06:29 

2012Nov06 

56237.61 

0.01 

S1S2E1E2W1W2 

02:38 

20130ct08 

56573.86 

0.66 

S1S2E1E2W1W2 

08:38 

2013Oct21 

56586.76 

0.76 

S1S2E1E2W1W2 

06:14 


vantage of the dense and complete ( u , v)-coverage from MIRC 
to create an image using aperture synthesis and to constrain the 
inclination angle of the disk. 

First, to extract the companion location, we fitted a bi¬ 
nary model to each of the ten best observing epochs with six- 
telescope MIRC data. The model for each epoch has 11 param¬ 
eters: fraction of H -band light from the primary, uniform disk 
diameter for the primary, fraction of light from the companion, 
uniform disk diameter for the companion, separation and posi¬ 
tion angle of the companion compared to the primary, Gaussian 
model for disk emission with three parameters describing the 



Fig. 3. Image reconstruction made with MIRC data in the H band. 
Here we have separately modeled the contributions from the Be star 
and the companion, and this image represents the emission from the 
disk alone, contributing about ~29% of the total H -band flux. We have 
marked the location of the companion at the ten epochs reported here 
with error ellipses and show the best-fit circular orbit (see text). The 
colored band shows the results from 100 bootstrap fits to illustrate the 
errors on our orbital solution. East is left and north is up. 


symmetric component (FWHM of the major axis, FWHM of the 
minor axis, position angle on the sky of the major axis), and two 
parameters to describe any asymmetric disk emission (skew an¬ 
gle and skew fraction; see description in Monnier et al. 2006; 
Che et al. 2012b). In our case, the uniform disk diameter (UDD) 
of the primary was fixed at 0.3 mas and the UDD of the com¬ 
panion was fixed as 0.056 mas, on the basis of spectrophotomet- 
ric estimations. We might expect the skewness to change from 
epoch to epoch (as was found for S Sco by Che et al. 2012b). 
This is indeed confirmed by our results at each epoch, but unfor¬ 
tunately, no conclusion could be drawn. Based on our fitting, we 
report the median and one-sigma range of the H -band disk pa¬ 
rameters in Table 5, while the locations of the companion for all 
epochs are collected in Table 6. The typical reduced^ 2 achieved 
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during fitting was ~1.9. A synthetic view of the MIRC data is 
presented in Fig. 2. 

We used the MACIM imaging algorithm (Ireland et al. 2006) 
along with the dark energy regularizer to carry out image recon¬ 
struction by combining all the MIRC data together. Since the 
primary star contains -70% of the total flux, it is important to re¬ 
move this component explicitly as part of the image reconstruc¬ 
tion, as explained in detail in Che et al. (2012b). We accounted 
for the moving secondary based on an orbital model, although 
this has only a small effect since the companion only contributes 
- 1.5% of the total light in the system. Figure 3 shows our image 
reconstruction of the disk alone, along with the locations of the 
secondary star for each epoch. The dynamic range appears to be 
noticeably smaller than most MIRC imaging (background noise 
is about ten times below peak emission), but this is because the 
central star has been removed from this image through model¬ 
ing. If we were to add back the central star, we would find a 
dynamic range of -100 in this image. The reduced^ 2 for this 
image is -2.0, only slightly worse than found using disk mod¬ 
eling for each individual epoch. We expect some time variation 
between epochs, and so it is not too surprising that the^ 2 is not 
as low as typically found when imaging with MIRC. 

We briefly mention some known limitations of our method 
and point out future work. First, the disk model is quite sim¬ 
plistic. As discussed elsewhere (e.g., Kraus et al. 2012; Che 
et al. 2012b), the gas disk is expected to have a hole in the center 
where the star is, and we also expect significant obscuration of 
the central star at this high-inclination angle. Furthermore, the 
imaging suggests a Gaussian disk profile, which is probably not 
a realistic model. The asymmetric component of the disk emis¬ 
sion changes the photo center of the primary object, thus adding 
errors to our estimate of the secondary separation vector. Indeed, 
it seems the companion location of 2012Nov06 is physically im¬ 
possible, and our estimate might have been affected by changing 
asymmetries in the inner disk that were not accurately modeled 
by our method. 

Table 5. Basic parameters of the H -band disk model from MIRC obser¬ 
vations. The component fluxes are indicated by / a , f d , and f b for the Be 
star, disk, and companion, respectively. 


Model Parameter 

Value 

fa 

0.70±0.06 

fd 

0.29±0.06 

fb 

0.015±0.003 

FWHM major (mas) 

1.44±0.41 

FWHM minor (mas) 

0.22±0.12 

PA of major axis (deg) 

-64±3 

Disk inclination (deg) 

82±4 


5.2. Orbital solution for the companion 

Based on the radial velocity measurements presented in Sect. 3 
and on the companion positions determined in Sect. 5.1, we 
performed a joint astrometric and radial velocity solution of 
the orbit. We adopted the period from Bozic et al. (1995), 
P = 126.6731 ± 0.0071 d, in the first combined fit. By combin¬ 
ing the epochs of the ascending node of the companion (equal to 
the velocity maximum for the secondary and the velocity min¬ 
imum for the Be star primary) from Bozic et al. (1995), Gies 
et al. (1998) (119 orbital cycles later), and our new solution (166 
orbital cycles later), we then derived improved estimates for the 
period and epoch. The revised period P = 126.6982 ± 0.0035 d 


Table 6. Position of the companion as a function of the RJD as detected 
in the various MIRC datasets. Separation is given in mas and position 
angle in degrees. Columns 5-7 give the uncertainty ellipse. 


RJD 

Phase 

Separation 

PA 

0~ma jor 

minor 

PA(cr) 

55807.91 

0.62 

4.50 

126.7 

0.07 

0.06 

84 

55832.86 

0.81 

2.56 

268.2 

0.09 

0.06 

18 

55852.87 

0.97 

5.75 

293.1 

0.11 

0.10 

12 

55853.76 

0.98 

5.83 

293.5 

0.14 

0.10 

298 

56158.02 

0.38 

4.34 

104.0 

0.10 

0.09 

330 

56185.81 

0.60 

4.90 

124.8 

0.14 

0.10 

318 

56231.77 

0.96 

5.73 

292.3 

0.10 

0.10 

291 

56237.61 

0.01 

6.08 

306.6 

0.18 

0.12 

292 

56573.86 

0.66 

3.39 

134.4 

0.26 

0.17 

340 

56586.76 

0.76 

1.27 

226.6 

0.28 

0.10 

309 



Fig. 4. Radial velocity measurements of the primary star along with 
predicted orbit from the joint MIRC+RV fit (see Table 7). Open sym¬ 
bols represent Her wing bisector measurements, while filled symbols 
represent HST measurements from Gies et al. (1998). 1 - cr error bars 
have been fixed to +8.6 km s -1 as explained in the text. The width 
of the solid line illustrates the errors on our orbital solution estimated 
from 100 bootstrap fits. The RV is negative at zero phase because we 
have chosen an astrometric orbital ephemeris of the secondary relative 
to the primary, which is opposite in sign to the radial velocity conven¬ 
tion, where the velocity of the primary star is measured relative to the 
center-of-mass. 

is 3 cr larger than that found by Bozic et al. (1995), but we use 
this revised period and epoch (r RV m i n = RJD 56110.00 ± 0.12) 
to reference the orbital phase in the rest of this paper. 

We fixed the period to P = 126.6982 and assumed a circular 
orbit for our final solution. To assess the errors in the orbital 
parameters, we carried out a bootstrap analysis, for which we 
resampled the astrometric and RV data to create 100 synthetic 
data sets. Initially, we found the reduced^ 2 for the RV data to 
be 69, clearly revealing systematic errors not accounted for in the 
original RV error analysis. In response to this, we treated all data 
with uniform error bars scaled to achieve a reduced x 2 of unity. 
The corresponding errors become thus ±8.6 km s -1 (resp. ±1.35 
km s -1 and ±3 km s -1 ) for the H a set (resp. the primary and 
secondary HST data set from Gies et al. 1998). As mentioned 
earlier, we also chose to remove the epoch of 2012 Nov 06 from 
our astrometric fit because it has a much larger discrepancy with 
the orbital prediction than any other observation, and it must be 
an outlier due to unmodeled disk substructure or a calibration 
problem (we have included it in Fig. 3 for reference). 
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D. Mourard et al.: Spectral and spatial imaging of the Be+sdO binary p Persei 


Table 7. Orbital elements of the p Persei binary system determined 
using interferometric and radial velocity measurements together, as¬ 
suming a circular orbit. Note our orbital convention for the ephemeris 
here is referencing the location of the secondary relative to the primary 
(thus the primary star attains its most negative radial velocity at epoch 

Try min)- 


Parameter 

Value 

Try min (RJD) 

56110.03+0.08 

P (d) 

126.6982 (fixed) 

a (mas) 

5.89+0.02 

e 

0 (fixed) 

m 

77.6±0.3 

con 

0 (fixed) 

no 

-64.3+0.3 

y (km s -1 ) 

-2.2+0.5 

K a (kms -1 ) 

10.2+1.0 

K b (kms -1 ) 

81.5+0.7 

M a+b (M 0 ) 

10.8+0.5 

M a (M 0 ) 

9.6+0.3 

Mb (M 0 ) 

1.2+0.2 

d(pc) 

186±3 


The result of our best-fit model along with 100 bootstraps is 
shown in Fig. 3, the RV curve is shown in Fig. 4. The results on 
the orbital parameters, masses, and distance appear in Table 7 
and are discussed in Sect. 7. 

6. Visible data and disk properties 

6.1. Visible broadband data and model fitting 

Using the VEGA wide spectral band measurements, we first 
performed a model fitting using the software LITpro 5 (Tallon- 
Bosc et al. 2008). We first used the three bands containing 
continuum and then, separately, the band containing the H a line 
as defined in Table 3. 

In a first approach, we just used a simple uniform disk model. 
We found a diameter in the continuum of 0.343 ± 0.006 mas and 
in the line of 0.423 ± 0.012 with a reduced^ 2 of 9.3 and 14.4, 
respectively (Table 8, first section). The larger equivalent diam¬ 
eter in the line indicates the presence of an extended structure, 
corresponding to the line emission associated with the disk. The 
residuals clearly indicate that the size of the object varies with 
position angle, and thus we used a second model with the addi¬ 
tion of an elongated Gaussian disk. The results are presented in 
Table 8 (second section). We tried to detect the companion in the 
VEGA data but, despite the very good ( u , v) coverage, the data 
quality at low visibility level does not permit this direct detec¬ 
tion. We therefore decided to adopt the position of the compan¬ 
ion found by MIRC, but consider its flux as a free parameter. To 
avoid the complication introduced by the change in the compan¬ 
ion’s position, we only used the data of 2011 October 18 and 19 
(RJD=55852 and 55853) for this analysis. 

The results clearly show the significant departure from a 
spherical object, both in the continuum and in the band con¬ 
taining H a. As evidenced by the evolution of the reduced x 1 , 
an improvement is found with a model consisting of a uniform 
disk + an elongated Gaussian disk (Table 8, third section), and 
the results are even better if we consider the companion (Ta¬ 
ble 8, fourth section). The difference in size and orientation of 

5 LITpro software available at http://www.jmmc.fr/litpro 


Table 8. p Persei model-fitting results for the continuum and the 
band containing Her. The different models are 1) UD=equivalent uni¬ 
form disk, 2) GD=elongated Gaussian disk, 3) UD+GD=central uni¬ 
form disk of flux FI + elongated Gaussian disk of flux F2, and 4) 
UD+GD+C=same as 3) + a companion of flux F3. 



Continuum 

Ha 

UD 




9.3+0.1 

14.4+0.1 

e 

0.343+0.006 0.423+0.012 

GD 




5.4+0.1 

7.2+0.1 

@min 

0.141+0.006 0.182+0.010 

elong 

2.34+0.14 

2.32+0.18 

PA 

-78+1 

-72+2 

UD+GD 



+ 

4.9+0.1 

4.7+0.1 

FI 

0.47+0.13 

0.45+0.08 

6 

0.1+0.1 

0.1+0.1 

F2 

0.53+0.13 

0.55+0.08 

@min 

0.19+0.02 

0.29+0.02 

elong 

2.64+0.49 

2.89+0.29 

PA 

-76+1 

-69+2 

UD+GD+C 



x 1 

4.3+0.1 

4.2+0.1 

FI 

0.44+0.06 

0.43+0.07 

e 

0.1+0.1 

0.1+0.1 

F2 

0.53+0.06 

0.54+0.09 

@min 

0.18+0.01 

0.28+0.02 

elong 

2.67+0.21 

2.92+0.30 

PA 

-76+1 

-69+2 

F3 

0.033±0.006 0.029±0.008 


the elongated object between the continuum and H a is the sig¬ 
nature of the difference between the disk emitting in the H a line 
and the continuum, and the elongated object potentially consists 
of a rotationally distorted photosphere plus the contribution from 
the disk. In fact, a large portion of the stellar photosphere is ob¬ 
scured by the disk, but this fraction is not the same in the con¬ 
tinuum and in the H a band. We also note that the position angle 
value found in the emission line region, where the disk is dom¬ 
inating the flux, is closer to the position angle found in the H 
band (PA = -64° ±3°) and is compatible within 2 cr despite the 
difference of wavelength. 

6.2. Spectrally resolved Ha data and model fitting 

Differential visibilities are estimated through the standard 
VEGA procedure (Mourard et al. 2009). For the purpose of the 
(p Persei work, we have established the absolute orientation of 
the differential phases in any of the CHARA+VEGA configu¬ 
rations (Mourard et al. 2012). Some examples of differential 
measurements are presented in Fig. 5, bottom part. The photon¬ 
counting detectors of VEGA suffer from a local saturation effect 
in presence of bright spectral lines that attenuates the line inten¬ 
sity. This nonlinearity in the photometric response of the detec¬ 
tor does not generate bias in the derived complex visibilities as 
explained by Delaa et al. (2011). 

To model the spectrally resolved visibility and differential 
phase in the H a emission line, we used the simple kinematic 
model developed for fast model-fitting of a geometrically and 
optically thin equatorial disk in rotation described in detail by 
Delaa et al. (2011) and Meilland et al. (2012). 
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The model parameters can be classified into four categories: 

1. The stellar disk parameters: stellar polar radius (R*), stellar 
flattening (f), distance (d), inclination angle (/), and major- 
axis position angle (PA). 

2. The kinematic parameters: rotational velocity (V rot ) at the 
disk inner radius (i.e., photosphere) and the rotation power- 
law index (J3 ). 

3. The disk continuum parameters: disk FWHM in the contin¬ 
uum (< a c ) and disk continuum flux normalized by the total 
continuum flux ( F c ). 

4. The disk emission line parameters: disk FWHM in the line 
(a\) and line equivalent width (EW). 


Table 9. Values of the VEGA best-fit kinematic model parameters. 


Parameter Value Remarks 


R* 

f 

d 

i 

PA 

Global geometric parameters 

7+2 R o 

1.4+0.3 

186 pc this work 

78°+10° 

-68°+ 5° 

V ro , 

Global kinematic parameters 

450 +50 km s -1 

P 

0.5+0.1 Keplerian rotation 

a c 

Continuum disk geometry 

3.4+1.7 D* = 1.2+0.6 mas 

F c 

0.14+0.2 

ana 

H a disk geometry 

5.9+1.7 D* = 2.1+0.6 mas 

EW Ha 

19.0+4.0 A 


The model produces a data cube consisting of narrow spec¬ 
tral band images through the Her line at the instrumental spec¬ 
tral resolution. Visibilities and differential phases for our base¬ 
line configurations are then extracted from the data cube and are 
compared to the observation using the standard^ 2 method. 

Except for the distance, the ten parameters are free and can 
vary within a range of realistic values. The parameter values 
for the best-fit model are presented in Table 9. The reduced^ 2 
for this model is 4.0. Interestingly, with this more representative 
model of the system geometry, the position angle and inclination 
of the disk is fully compatible with the one found with the MIRC 
imaging. 

The 0.65yum continuum-squared visibilities for the best-fit 
model are plotted as a function of the spatial frequency in Fig 5, 
top part. Some examples of differential visibilities and phases 
for the best-fit model for various baselines between 32 and 278m 
are plotted in Fig 5, bottom part. As already explained in De- 
laa et al. (2011), phase variations for baselines overresolving the 
object (i.e., 60m in the major-axis orientation and 200m in the 
minor-axis one in the case of cp Persei) are only poorly fitted. 
The jumps seen for some baselines might be either due to inho¬ 
mogeneity in the disk or to VEGA data reduction biases at very 
low visibility. We also note that some of the very short baselines 
show an increase of visibility instead of a drop in the H a line, 
which corresponds to a smaller object in the line. This effect 
might be due to the high opacity of the disk seen edge-on and 
can only be modeled using a radiative transfer code. 



1.5 


B=32.7m PA=2.8° 

B = 65.8m PA=42.5° 

B=98.4m PA= 132.4° 

B= 173.7m PA=172.3° 





B=206.1m PA= —176.3° 

r ; 

B = 217.4m PA=30.6° 

; 

f ; 

B=241.0m PA=47.1° 

: V r 

B=278.5m PA= 167.7° 


\ 

\ 


650 655 660 665 

A (nm) 

650 655 660 665 

A (nm) 

650 655 660 665 

A (nm) 

650 655 660 665 

A (nm) 


Fig. 5. Top: 0.65//m continuum visibility for the best-fit model is 
plotted as a function of the spatial frequency, whereas the blue zone 
represents the extent of the model visibility for the various orientations 
of the baseline. Bottom: examples of differential visibilities and phases 
(in red + error bar in yellow) for the best-fit model (in blue) for various 
baselines between 32 and 278m. 


6.3. Spectrally resolved image reconstruction around Ha 

We used the MiRA software (Thiebaut & Giovannelli 2010) 
combined with an improved self-calibration method (Millour 
et al. 2011) to reconstruct wavelength-dependent images of <p 
Persei in the visible wavelength range. The data set contains the 
broadband V 2 defined in Table 3 and the wavelength-dependent 
differential visibilities and phases, but no closure phases, as pre¬ 
viously said. Therefore, the self-calibration method had to be 
modified. The different steps of our method, and the additions 
made in this work, are presented in Appendix A. The so-called 
SELF-CAL software and its documentation are available online 6 . 

The image reconstruction leads to an (x, y. A) image cube that 
covers a field of view of 6.4 mas with 64x64 pixels, a wavelength 
range between 646 nm and 666 nm with a spectral resolution of 
« 1640. This translates into ^10 spectral channels in the H a 
emission line. 

The images in Fig. 6 exhibit little variation of morphology 
throughout the continuum. They display three prominent fea¬ 
tures: 


6 Available at http://self-cal.oca.eu 
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D. Mourard et al.: Spectral and spatial imaging of the Be+sdO binary ip Persei 
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Fig. 6. Top row: sequence of intensity maps (6.4x6.4mas) as a function of wavelength with north up and east left as imaged by VEGA (upper 
sequence), as computed from the kinematic model in Sect. 6.2 (middle sequence), and as imaged through a simulated observation of the above 
model (lower sequence). Middle row: at the left, we present the continuum image convolved to the diffraction limit of VEGA (0.2 mas) together 
with contours at 50, 5, 0.5 and 0.05% of the peak intensity to show the disk contour. The middle and right panels show the selected zones to 
extract the spectra presented below together with an unconvolved continuum image of the system. Bottom row: extracted spectra for the different 
spatial components: 1) Bottom middle: extracted fluxes (top: total flux, pink absorption feature: “star” zone, blue emission feature: “disk” zone, 
and dashed spectrum: all other contributions), 2) bottom right: decomposition of the “disk” zone: the top spectrum is the total flux. The blue and 
red lines correspond to the blue- and redshifted emission features, while the central zeroshifted emission feature comes from all other parts of the 
“disk” zone. 


- an elongated central source (the “star”) that can be contained 
in an ellipse with « 0.6 x 0.4 mas major and minor axes, and 
« 20° position angle for the minor axis, 

- an extended flattened source (the “disk”) whose extension is 
weakly constrained. We can set an upper limit to its size to 
an ellipse « 3.9 x 1.8 mas, and « 20° position angle of the 
minor axis, 

- a significant amount of unconstrained flux (« 15% of the 
continuum flux, throughout all wavelengths), which could 
be linked to the companion star or to image-reconstruction 
artifacts. 

The images in the emission line are very different from 
those in the continuum, showing emission that sweeps across 


the “disk” feature. In the upper part of Fig. 6, we compare the 
spectral images and the maps deduced from the model presented 
in the previous section. We also provide a set of reconstructed 
images from simulated data based on the kinematic model. 

To extract information from these images from milli- 
arcsecond integral-field spectroscopy, we multiplied the im¬ 
ages by the system spectrum (as MiRA provides total-flux- 
normalized images) and selected different zones with binary 
masks to extract their spectra. The result of this process is shown 
in the bottom part of Fig. 6. 

The “star” zone exhibits an absorption spectrum, character¬ 
istic of a stellar spectrum. The “disk” zone bears most of the 
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emission of the system, while the unconstrained flux has a flat 
spectrum. 

The results of the VEGA imaging and model fitting are dis¬ 
cussed in general in Sect. 7. 

7. Discussion 

7.1. Fundamental stellar parameters 

We first consider the dimensions of the binary system. The as¬ 
trometric measurements from MIRC (Table 6) were combined 
with the radial velocity measurements for the two components 
(Table 2) in a joint solution for the binary orbital elements and 
distance that is summarized in Table 7. Our derived values of the 
stellar masses, M a = 9.6 ± 0.3M© and M b = 1.2 + 0.2A/©, agree 
well with those determined previously by Gies et al. (1998). Our 
distance estimate d = 186 ± 3 pc is about 4.2 cr closer than the 
estimate from the revised Hipparcos distance or 1. lcr closer than 
the original Hipparcos distance. We note, however, that p Persei 
is flagged as ’’unsolved variable” in the Hipparcos catalog; we 
could then expect an underestimated error for the parallax in the 
revised Hipparcos catalog. 

A knowledge of the distance to p Persei has consequences 
for the size of the Be star. We can estimate the stellar radius 
from the predicted angular size and distance. We first consider 
the angular size from a comparison from a of the observed and 
predicted photospheric fluxes, a task that is best accomplished in 
the ultraviolet part of the spectrum where the circumstellar disk 
flux is negligible. Touhami et al. (2013) fit the observed UV 
spectrum of p Persei using a Kurucz flux model for the adopted 
temperature (from Fremat et al. (2005)) of the Be star and a stan¬ 
dard extinction law (ignoring the small flux contribution of the 
hot companion), and they derived a limb-darkened, angular di¬ 
ameter for the Be star of 0.235 ± 0.008 mas. This is similar to the 
minor-axis diameter derived from the VEGA broadband obser¬ 
vations (0.18 + 0.01 mas, Table 8). Then, for an adopted distance 
of 186 pc, the stellar radius is estimated to be R a /Ro = 4.7 ± 0.4. 
This is somewhat lower but consistent within error bars to that 
given in Table 9 and prorated for the adopted lower distance 
( R a /RQ = 5.9 ± 1.7). The difference may be partially due to the 
neglect of UV flux scattered away from the line of sight by elec¬ 
tron scattering in the inner disk so that the stellar UV flux (and 
hence angular size) is somewhat higher than estimated from the 
observed UV spectrum. Furthermore, the estimate of stellar di¬ 
ameter from the VEGA measurements may be influenced by the 
assumed flux distribution in the disk. 

The mass and radius of the Be star match the mass - ra¬ 
dius relationship for main-sequence stars derived from eclips¬ 
ing binaries (Torres et al. 2010). Furthermore, the effective 
temperature T eff = 29.3 kK (Gies et al. 1998) and luminos¬ 
ity lo gL a /L Q = 4.16 ± 0.10 (for R a /Ro = 4.7) are consis¬ 
tent with those of eclipsing binary stars of similar mass (Tor¬ 
res et al. 2010). The Be star’s rotation rate, however, is excep¬ 
tional. Adopting the estimate of projected rotational velocity 
from Fremat et al. (2005), V sin i = 462 ± 33 km s -1 , we derive 
a ratio of equatorial velocity to orbital velocity at the equator 
(Rivinius et al. 2013) of (V sin i/ sin 0/Vcrit = 0-93 ± 0.08 for the 
mass and radius given above and an inclination from Table 5 (as¬ 
suming coalignment of the spin and orbital vectors; see below). 
Note that we have tacitly assumed that the radius derived from 
the UV flux and the distance corresponds to the star’s polar ra¬ 
dius, whereas it may represent some value between the polar and 
(larger) equatorial radius. For example, if our radius estimate 
corresponds to the equatorial radius in a near critical rotation 


star, then R ±po iar = R a / 1.5 and (V sin i/ sin 0/Vcrit = 0-26 ± 0.08. 
On the other hand, the measured projected velocity may under¬ 
estimate the actual value of V sin i because of gravity darkening 
(Townsend et al. 2004). Despite these uncertainties, the clear 
conclusion is that the Be star is rotating very rapidly. The com¬ 
panion’s rotational velocity is low, v rot sin/ < 10 km s -1 (Gies 
etal. 1998). 

Gies et al. (1998) derived a radius ratio of Rb/R a = 0-20 
from the observed flux ratio in the ultraviolet spectrum. Adopt¬ 
ing the Be star radius given above, the radius and luminosity of 
the hot companion are Rb/Ro = 0.94 ± 0.09 and logZ^/L© = 
3.80 ± 0.13 (for T eff = 53 + 3 kK; Gies et al. (1998)). Thus, 
although the Be star is five times larger than the sdO compan¬ 
ion, it is only twice as luminous. In the Rayleigh-Jeans tail of 
the spectrum in the near-IR, we would expect that the surface 
flux varies with effective temperature, so the flux ratio would 
be Fb / F a « ( Tb/Ta)(Rb/R a ) 2 - Based upon the temperature and 
radius ratios estimated by Gies et al. (1998), the predicted opti- 
cal/IR flux ratio is 8%. This value is similar to the observed flux 
ratio in the optical (F3/F1 = 7.5 + 1.7%; Table 8), but higher 
than found in the H band (/ 2//1 = 2.1 + 0.5%; Table 5). We 
speculate that the discrepancy is due to obscuration of the sdO 
companion by circumstellar gas that becomes more significant 
at longer wavelength (see below). 

7.2. Birthplace and age 

We next consider the possible birthplace and age of p Persei. Our 
derived distance is consistent with that of the nearby a Per clus¬ 
ter (Melotte 20; Per OB3) at 111 pc (de Zeeuw et al. 1999), and p 
Persei is located at a similar angular separation (15.9°) from the 
cluster center as some of the a Per “halo stars” identified by de 
Zeeuw et al. (1999). Furthermore, the proper motion of p Persei 
( p a cos 6 = 24.6 mas yr -1 , p 6 - -14.0 mas yr -1 ; van Feeuwen 
2007) is similar to that of the a Per cluster (p a cos S = 22.7 
mas yr -1 , p 5 = -26.5 mas yr -1 ; Fodieu et al. 2012). Fikewise, 
the systemic radial velocity of p Persei (-2.2 km s -1 ; Table 7) is 
similar to the mean value for the a Per cluster (-1.0 km s -1 ; de 
Zeeuw et al. 1999). Thus, the distance, position, and kinematical 
properties of p Persei suggest that it may be an outlying member 
of the a Per cluster. If so, then the age of p Persei is the cluster 
age, « 52 Myr (Kharchenko et al. 2005; Makarov 2006; Fyu- 
bimkov et al. 2010). This is much older than the main-sequence 
lifetime of « 25 Myr for a 10M© mass star like the Be primary 
(Brott et al. 2011), which suggests that the Be star attained its 
high mass relatively recently. 

The age of the a Per cluster is consistent with the binary 
evolutionary scenario for p Persei. The p Persei binary probably 
began its life as a pair of 6M© + 5 M© stars with a shorter or¬ 
bital period in the range of 7 d (Vanbeveren et al. 1998a) to 14 d 
(Vanbeveren et al. 1998b). As the initially more massive star 
completed core-hydrogen burning, it expanded and began mass 
and angular momentum transfer to the companion. The binary 
separation increased following mass ratio reversal, reaching its 
current state of wide separation between the stripped-down rem¬ 
nant of the mass donor (the hot subdwarf star) and the rapidly 
rotating mass gainer (the Be star). If we assume that the mass 
transfer was approximately conservative (i.e., negligible mass 
lost from the system), then the progenitor of the hot subdwarf 
had a mass at least as high as half the current total mass of the 
system (5.4 M©; Tab. 7). Indeed, given the current advanced evo¬ 
lutionary state of the donor as a hot subdwarf star, it must have 
had an original mass in excess of that for any H-buming star still 
in the a Per cluster. The most massive star in the a Per clus- 
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ter (aside from ip Persei) is a Per itself, with a mass estimated 
in the range of 6.7M© (Makarov 2006) to 7.3M© (Lyubimkov 
et al. 2010). Thus, the original mass of the donor star was prob¬ 
ably Mb orig ~ 7Af©, so that the mass gainer (Be star) accreted 
most of the difference Mb or ig - M b ~ 6M© during mass transfer. 
This wholesale transformation of the mass gainer led to its cur¬ 
rent high mass and fast spin, and the Be disk found today repre¬ 
sents the Be star’s means to shed its excess angular momentum. 
We now know of two other Be plus sdO binaries, FY CMa and 
59 Cygni, that share many of the properties of ip Persei (Peters 
et al. 2013), and it is probable that other hot companions of Be 
stars await discovery (Koubsky et al. 2012, 2014). 

7.3. Properties of the Be-star disk 

We next consider the properties of the Be star disk. Interferomet¬ 
ric observations have now resolved the disk in the continuum at 
several wavelengths. The angular size derived depends in part on 
the adopted model for the flux distribution in the sky, but we list 
in Table 1 the results for an assumed uniform disk star and Gaus¬ 
sian elliptical disk. The Gaussian FWHM along the disk major 
axis 6 ma j appears to increase with wavelength from R through 
H and K' bands as expected because of the increase in hydro¬ 
gen free-free opacity with wavelength that causes the optical 
depth unity radius to appear larger when moving into the near- 
infrared (Touhami et al. 2011). The H a diameter measured from 
the VEGA observations is somewhat smaller than that derived 
from earlier H a measurements and is comparable to or lower 
than the contemporaneous near-IR continuum estimates. This is 
surprising given that the H a sizes are generally larger than the 
near-IR continuum sizes in most Be stars (Touhami et al. 2013). 
We note, however, that the H a emission equivalent width at the 
time of the VEGA observations (W,\ = -25 A) was much weaker 
than at some times in the past ( W\ = -46A in 2001), and if 
the disk size scales with the emission strength (Grundstrom & 
Gies 2006), then the disk dimensions may have been larger at 
earlier epochs. It is possible that gravitational interactions with 
the companion limit the physical boundary of the gaseous disk, 
as occurs in many Be X-ray binaries (Negueruela & Okazaki 
2001 ). 

The interferometric observations indicate that the spin angu¬ 
lar momentum axis of the disk is coaligned with the orbital an¬ 
gular momentum axis. The position angle of the long axis of the 
projected disk (PA in Tables 1, 5, 8, and 9) agrees with the lon¬ 
gitude of the ascending node in the astrometric orbital solution 
(Q in Table 7). The position angle of the disk normal (PA + 90°) 
agrees with the intrinsic polarization angle, as expected by Quir- 
renbach et al. (1997). Furthermore, the disk inclination derived 
from the ratio of projected minor to major axis (i = 75° ±9° from 
r = cos i in Tables 1 and 9) is consistent within errors with the 
astrometric orbital inclination (i = 77?6±0?3 in Table 7). Finally, 
the sense of the disk rotation (moving away in the northwestern 
part; Fig. 6) is the same as the sense of orbital motion (Fig. 3). 
Thus, we find that the disk and orbital angular momentum vec¬ 
tors are coaligned, verifying the prediction from the binary in¬ 
teraction model that the spin-up of the Be star was caused by 
mass and angular momentum transfer from the donor compan¬ 
ion. Note that in general we expect that the inclination derived 
from the projected shape of the disk will be smaller than actual 
because of disk flaring that increases the apparent size of the 
minor axis (Grundstrom & Gies 2006), so the small difference 
we find between inclination estimates from the disk and orbit is 
probably consistent with some disk flaring. The increase in disk 


thickness with radius (flaring) means that there will be a range 
in inclination below 90° where the disk may occult our view of 
the Be star photosphere, and this may explain the presence of 
narrow “shell” lines that are observed in the spectra of some Be 
stars (Rivinius et al. 2013). ip Persei has a well-observed shell 
line spectrum, so our inclination estimate would appear to indi¬ 
cate that the disk half-opening angle may amount to as much as 
90° - i = 12° in this case. Curiously, Hanuschik (1996) used a 
statistical study of the numbers of Be-shell and normal Be stars 
to arrive at a similar half-opening angle of 13°. The agreement 
may be coincidental, and detailed models of disk line formation 
show that shell line spectra may form over a wider range below 
i = 90° (Silaj et al. 2014). 

Next, we offer a few comments about the surface brightness 
distribution of the Be disk. Similar to other investigators, we 
have used Gaussian elliptical models for surface brightness be¬ 
cause they generally provide good fits of the interferometric vis¬ 
ibility. However, models of the spatial appearance of disk flux 
often predict more complicated distributions (Fig. 6), and in par¬ 
ticular, in models where the circumstellar gas becomes optically 
thick, the surface brightness distribution appears more or less 
constant in the parts of the disk close to the photosphere (Gies 
et al. 2007). Thus, for example, the relatively small minor axis fit 
to the MIRC observations (Table 5) may be partially due to the 
non-Gaussian appearance of the foreshortened disk brightness 
distribution (Fig. 3). It is probably unrealistic to assume that the 
disk is azimuthally symmetric, because many Be stars including 
ip Persei display long-lived asymmetries in their emission line 
profiles that probably form in one-armed spiral structures (Riv¬ 
inius et al. 2013). The H a profile observed over the time span 
of the CHARA observations displayed a stronger red than blue 
peak (see Fig. 8 of Silaj et al. (2010)), which indicates that some 
asymmetry in the H a disk surface brightness did exist at this 
time. We note that there is a hint in the H a image reconstruc¬ 
tions shown in Fig. 6 that the northwestern (redshifted) exten¬ 
sion of the disk is slightly brighter than the southeastern exten¬ 
sion (blueshifted). The VEGA observations were obtained near 
orbital phase 0.0 when the companion appears near its north¬ 
western, maximum elongation, and the part of the disk facing 
the companion may be brighter (Hummel & Stefl 2001). Thus, 
the H a interferometry offers tentative evidence for the presence 
of the surface brightness asymmetries expected from disk den¬ 
sity structures and heated zones. 

Finally we note that reverse mass transfer from the Be disk 
to the hot subdwarf may occur. The fact that the observed cir¬ 
cumstellar disk radius appears smaller than the semimajor axis 
(Fig. 3) does not necessarily imply that disk gas never reaches 
the vicinity of the companion. The interferometric observations 
record the spatial distribution of the brightest parts of the disk 
emission, and it is possible that lower density and fainter disk 
gas does extend to the orbit of the companion as it slowly moves 
outward from the Be star. Indeed, Hummel & Stefl (2003) ar¬ 
gued that the He I 6678A emission extends to the outer edge 
of the disk at RR 0C he/^ = 0-56. We would expect that some of 
this gas will be captured by the gravitational pull of the com¬ 
panion to form an accretion disk around the subdwarf star. The 
double-peaked appearance of the He II 4686A emission that 
forms near the companion may be evidence of such an accretion 
disk (Poeckert 1981). It is possible that gas around the compan¬ 
ion may obscure our view of the hot subdwarf star itself, and 
the attenuation would increase with wavelength if free-free pro¬ 
cesses dominate the continuum opacity. This may explain why 
the companion appears fainter than expected in the near-IR. 
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8. Conclusion 

Thanks to advanced visible and infrared interferometric instru¬ 
mentation at the CHARA Array, we have been able to obtain the 
first high angular resolution images of the Be binary p Persei, in¬ 
cluding the first spectrally resolved images in the visible (made 
with VEGA in four-telescope mode). These images show the 
highly inclined circumstellar disk surrounding the Be star, and 
the helium star remnant of the mass donor star is detected for the 
first time through MIRC observations of its orbital motion. 

The relative angular orbital motion of the hot companion was 
combined with new orbital radial velocity measurements of both 
components to derive a full, three-dimensional orbit of the binary 
system. We obtain masses of 9.6+0.3M© and 1.2+0.2M© for the 
Be star primary and subdwarf secondary, respectively. The dis¬ 
tance derived from the joint astrometric and spectroscopic solu¬ 
tion is 186+3 pc. We note that p Persei shares the same distance, 
kinematical, and evolutionary properties as members of the a Per 
star cluster, and if this association is confirmed, then the cluster 
age will place significant constraints on the initial masses and 
evolutionary mass transfer processes that transformed the p Per¬ 
sei binary system. The Be star (mass gainer) has a mass, ra¬ 
dius, and effective temperature that are typical of an early B-type 
main-sequence star, but it is rotating very rapidly (93 ± 8% of 
the critical rate) as a consequence of angular momentum gained 
during past mass transfer. The sdO secondary star (remnant of 
the mass donor) is fainter in the near-IR than predicted from the 
FUV flux ratio (Gies et al. 1998), due perhaps to obscuration by 
circumstellar gas in its vicinity. 

The circumstellar gas disk of the Be star primary was re¬ 
solved in the optical and near-IR H- band, and we estimated the 
disk flux contributions in each of the H a emission lines and the 
optical and near-IR continua. The apparent ellipticity and po¬ 
sition angle of the disk in the sky is similar in each band. The 
size of the disk in the continuum flux is larger at longer wave¬ 
length, presumably due to the increasing free-free optical depth 
of the disk at longer wavelength. The disk appears largest in the 
spectrally resolved H a observations (Fig. 6), but the visible flux 
appears to be confined well within the Roche lobe of the Be star. 
We find that the angular momentum vector of the disk is aligned 
with and shows the same sense of rotation as the binary orbital 
angular momentum vector. This agreement confirms the evolu¬ 
tionary prediction that the Be star was spun up by mass transfer 
in the binary and that this excess angular momentum is now be¬ 
ing shed into the Be star disk. 

The appearance of the Be disk is presented in the first non- 
parametric images reconstructed from the MIRC H -band obser¬ 
vations (Fig. 3) and from the spectrally dispersed VEGA obser¬ 
vations of H a (Fig. 6). These images reveal the brighter, inner 
parts of the disk, and there is some evidence in the spectrally dis¬ 
persed H a maps from VEGA of an asymmetry in the azimuthal 
light distribution of the circumstellar disk. The northwestern 
part of the disk that is moving away from us appears slightly 
brighter than the southeastern part that is approaching (Fig. 6), 
consistent with the relatively stronger red peak of the H a profile 
observed at this time. Such asymmetries probably originate in 
long-lived disk structures (perhaps a one-armed spiral instabil¬ 
ity), and long-baseline interferometric observations now appear 
to offer the means of exploring the time evolution of these disk 
processes. 
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Table 2. Radial velocity measurements 


Component 

RJD 

Vr (km s A ) 

Source and Observer 

Be Primary 

50023.5987 

-9.8 ± 1.5 

HST Gies etal. (1998) 

Be Primary 

50084.5151 

8.7 ± 1.5 

HST Gies etal. (1998) 

Be Primary 

50276.5235 

-12.8 ± 1.5 

HST Gies etal. (1998) 

Be Primary 

50327.3240 

4.1 ± 1.5 

HST Gies etal. (1998) 

Be Primary 

50371.3607 

-0.2 + 1.5 

HST Gies et al. (1998) 

Be Primary 

50371.4089 

-0.2 + 1.5 

HST Gies etal. (1998) 

Be Primary 

50709.6087 

1.9 + 1.0 

BeSS Neiner 

Be Primary 

50709.6193 

2.0+ 1.0 

BeSS Neiner 

Be Primary 

51056.9586 

-22.0+ 1.1 

KPNO Gies 

Be Primary 

51221.3465 

-0.2 + 0.8 

INAF Tomasella et al. (2010) 

Be Primary 

51419.9373 

-11.1 + 3.8 

KPNO Gies 

Be Primary 

51419.9390 

-11.3 + 3.5 

KPNO Gies 

Be Primary 

51420.9717 

-11.6 ±5.0 

KPNO Gies 

Be Primary 

51421.9875 

-14.5 + 2.6 

KPNO Gies 

Be Primary 

51421.9905 

-13.1 +2.7 

KPNO Gies 

Be Primary 

51430.9768 

-14.7 + 4.0 

KPNO Gies 

Be Primary 

51430.9805 

-14.8 + 4.2 

KPNO Gies 

Be Primary 

52237.3904 

-5.2 + 0.7 

BeSS Buil 

Be Primary 

52241.3562 

-10.1+0.7 

BeSS Buil 

Be Primary 

52253.3768 

-4.9 + 0.7 

BeSS Buil 

Be Primary 

52263.2743 

-2.3 + 0.4 

BeSS Buil 

Be Primary 

52284.3159 

-18.4 + 0.7 

BeSS Buil 

Be Primary 

52293.3072 

-18.0 + 0.6 

BeSS Buil 

Be Primary 

52303.2874 

-29.4 + 0.6 

BeSS Buil 

Be Primary 

52521.4100 

10.7 + 0.9 

BeSS Buil 

Be Primary 

52854.5309 

13.4 + 1.0 

BeSS Buil 

Be Primary 

52855.5299 

10.3 + 1.5 

BeSS Buil 

Be Primary 

53290.8068 

-0.7 + 1.2 

KPNO Grundstrom 

Be Primary 

53290.8076 

-0.6 + 1.4 

KPNO Grundstrom 

Be Primary 

53290.8087 

-0.0 ± 1.6 

KPNO Grundstrom 

Be Primary 

53291.7942 

-0.8 + 1.1 

KPNO Grundstrom 

Be Primary 

53291.7951 

-0.3 + 1.3 

KPNO Grundstrom 

Be Primary 

53292.8129 

-1.6 + 1.0 

KPNO Grundstrom 

Be Primary 

53292.8140 

-2.1 + 1.1 

KPNO Grundstrom 

Be Primary 

53294.7812 

-0.1 +2.2 

KPNO Grundstrom 

Be Primary 

53294.7851 

-0.1 + 1.7 

KPNO Grundstrom 

Be Primary 

54019.8264 

10.2 + 1.2 

KPNO Grundstrom 

Be Primary 

54020.8271 

8.1 +0.9 

KPNO Grundstrom 

Be Primary 

54021.8123 

6.6+ 1.2 

KPNO Grundstrom 

Be Primary 

54022.7978 

8.3 + 1.4 

KPNO Grundstrom 

Be Primary 

54023.8106 

7.0 + 0.9 

KPNO Grundstrom 

Be Primary 

54331.5338 

-24.6 + 2.1 

BeSS Desnoux 

Be Primary 

54358.5030 

-18.2 + 0.9 

BeSS Thizy 

Be Primary 

54393.5289 

5.6+ 1.2 

BeSS Thizy 

Be Primary 

54513.3741 

12.7 + 2.4 

BeSS Guarro 

Be Primary 

54672.6520 

3.2+ 1.0 

BeSS Buil 

Be Primary 

54761.7315 

-0.6 + 3.2 

KPNO Grundstrom 

Be Primary 

54793.4824 

8.3 + 1.1 

BeSS Guarro 

Be Primary 

54866.2563 

-27.1+2.1 

BeSS Guarro 

Be Primary 

55059.5401 

-8.5 + 3.5 

BeSS Desnoux 

Be Primary 

55060.5620 

-3.6 + 1.2 

BeSS Buil 

Be Primary 

55067.5022 

-9.5 + 1.4 

BeSS Thizy 

Be Primary 

55072.3554 

11.0 + 4.3 

BeSS Terry 

Be Primary 

55109.4944 

-12.6 + 2.8 

BeSS Buil 

Be Primary 

55117.3864 

-10.7 + 10.7 

BeSS Garrel 

Be Primary 

55119.3849 

-15.7 + 11.4 

BeSS Garrel 

Be Primary 

55203.3741 

-17.7 + 17.1 

BeSS Garrel 

Be Primary 

55398.4875 

11.1 + 1.3 

BeSS Terry 

Be Primary 

55417.5052 

8.0+ 2.5 

BeSS Desnoux 

Be Primary 

55455.4248 

-16.4+ 1.2 

BeSS Buil 

Be Primary 

55534.2922 

-1.9 + 1.6 

BeSS Guarro 
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Be Primary 

55559.4242 

-0.3 ± 2.0 

BeSS Garrel 

Be Primary 

55774.5908 

3.2 ±2.2 

BeSS Ubaud 

Be Primary 

55777.5094 

-2.6 ± 0.8 

BeSS Dubreuil 

Be Primary 

55785.5531 

2.5 ± 0.7 

BeSS Terry 

Be Primary 

55828.3888 

-17.2 + 0.6 

BeSS Terry 

Be Primary 

55831.5150 

-20.7+ 1.3 

BeSS Buil 

Be Primary 

55857.4198 

-31.3 + 0.6 

BeSS Buil 

Be Primary 

55879.6551 

-9.6 + 0.6 

BeSS Graham 

Be Primary 

55907.6115 

8.0+ 0.6 

BeSS Graham 

Be Primary 

55956.3393 

-12.3 + 2.1 

BeSS Garrel 

Be Primary 

56143.6007 

5.8 + 1.3 

BeSS Buil 

Be Primary 

56150.4413 

7.5 + 2.8 

BeSS Terry 

Be Primary 

56163.4058 

17.4 + 2.6 

BeSS Garrel 

Be Primary 

56167.4575 

18.1 + 1.8 

BeSS Terry 

Be Primary 

56204.5901 

0.8 + 2.2 

BeSS Favaro 

Be Primary 

56206.4024 

-0.4 + 0.7 

BeSS Buil 

Be Primary 

56260.5174 

-4.8 + 0.5 

BeSS Graham 

Be Primary 

56267.3228 

2.2+ 1.2 

BeSS Garrel 

Be Primary 

56270.4680 

4.7 + 2.0 

BeSS Guarro 

Be Primary 

56291.3623 

20.6 + 0.5 

BeSS Ubaud 

Be Primary 

56301.3892 

19.6 + 3.4 

BeSS Guarro 

Be Primary 

56490.5331 

-4.6 + 1.3 

BeSS Ubaud 

Be Primary 

56628.5588 

4.6+ 1.6 

BeSS Graham 

Be Primary 

56665.6322 

28.1 + 1.0 

BeSS Sawicki 

sdO Secondary 

43763.4936 

-114.3 + 3.5 

IUE 

sdO Secondary 

43837.9529 

62.5 + 3.6 

IUE 

sdO Secondary 

43839.8268 

48.4 + 2.6 

IUE 

sdO Secondary 

43844.0271 

35.0 + 3.2 

IUE 

sdO Secondary 

44151.5341 

-66.7 + 4.3 

IUE 

sdO Secondary 

44167.3172 

-8.6+ 3.8 

IUE 

sdO Secondary 

44174.0036 

7.6+ 3.8 

IUE 

sdO Secondary 

44328.2958 

81.8 + 4.5 

IUE 

sdO Secondary 

44332.4182 

86.0 + 4.0 

IUE 

sdO Secondary 

44485.0176 

-13.3 + 2.7 

IUE 

sdO Secondary 

44486.1504 

-16.9 + 2.7 

IUE 

sdO Secondary 

44511.1058 

-96.3 + 2.6 

IUE 

sdO Secondary 

44661.1698 

-82.7 + 2.6 

IUE 

sdO Secondary 

44824.9353 

66.5 + 2.6 

IUE 

sdO Secondary 

45237.9853 

28.6 + 3.2 

IUE 

sdO Secondary 

45238.0244 

26.5 + 4.4 

IUE 

sdO Secondary 

50023.5987 

78.2 + 1.1 

HST Gies etal. (1998) 

sdO Secondary 

50084.5151 

-77.6+ 1.1 

HST Gies etal. (1998) 

sdO Secondary 

50276.5235 

75.5 + 1.1 

HST Gies etal. (1998) 

sdO Secondary 

50327.3240 

-54.1 + 1.1 

HST Gies et al. (1998) 

sdO Secondary 

50371.3607 

-24.2+ 1.1 

HST Gies etal. (1998) 

sdO Secondary 

50371.4089 

-23.0+ 1.1 

HST Gies et al. (1998) 
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Appendix A: Visible-image reconstruction process 

Appendix A. 1: Self-calibration 

We used the software MiRA (Thiebaut & Giovannelli 2010) 
combined with an improved self-calibration method (Millour 
et al. 2011) to reconstruct wavelength-dependent images of y 
Persei in the visible wavelength range. The data set contains only 
broadband V 2 defined in Table 3 and the wavelength-dependent 
differential visibilities and phases, but no closure phases. There¬ 
fore, the self-calibration method had to be updated: the first 
step of image-reconstruction is made using the broadband V 2 
only. The first tentative map is therefore the same for all wave¬ 
lengths. Then we had to self-calibrate the visibility gains in ad¬ 
dition to the differential phases offsets. We outline below the dif¬ 
ferent steps of our method and the additions made in this work. 
We also refer to the review of Pearson & Readhead (1984) for 
an overview of the self-calibration methods in radio-astronomy. 
The so-called SELF-CAL software and its documentation are 
available online 7 . 


Pixel size The pixel size was selected thanks to 

(Z)pixel — A m \ n /£max/4. (A.4) 

The factor 4 comes from a double factor 2: the diffraction 
limit of the interferometer as defined as in Monnier et al. (2007), 
and an additional "super-resolution" factor 2. We have a max¬ 
imum baseline of 312 m, which translates into a 0.1 mas pixel 
size. 

In addition, for computing time considerations, we had to 
select a power of 2 for the pixel number. Therefore, we selected 
images with 64x64 pixels of 0.1 mas. 

Hyperparameter selection For the images presented here, we 
used the total variation regularization as suggested in Renard 
et al. (2011), a number of iterations per wavelength of 50, and 
we selected the hyperparameter value p = 10 5 thanks to the same 
method as described in Domiciano de Souza et al. (2014). 


1. Initial image reconstruction : we call /o(v,y, do) the first 
tentative image reconstruction using the V 2 (u, v, To) only, 
Iq(u,v,Aq) its Fourier transform, and dV(u, v, A), df(u,v,A) 
the differential visibilities and phases provided by VEGA. 

2. Self-calibration : in the self-calibration process we search for 
complex gain values g(u , v) that are a function of the (w, v) 
coordinates and have the following equalities: 

||/ 0 (w,v,d)|| = \\g(u,v)\\xdV(u,v,A) (A.l) 

arg[/ 0 (w, v, d)] = arg[g(w, v)] x dfi(u, v, d) (A.2) 


3. (x,y 9 A) Image reconstruction : these gain values are then ap¬ 
plied to the differential phases and visibilities, which become 
amplitude and phase of a complex visibility that is stored to 
a new OIFITS file (Pauls et al. 2005). This OIFITS file is fed 
again to the MiRA software to produce a (v, y, d) image. 

This new image serves as a basis to a new step of self¬ 
calibration. The self-calibration process can be repeated until 
a reasonable convergence is achieved. Here, we decided to arbi¬ 
trarily stop at the second step, as in Millour et al. (2011). 


Appendix A.2: Selecting proper image reconstruction 
parameters 

For the image reconstruction process, we had to select the field 
of view, the pixel size and the reconstruction hyperparameter p 
values. 


Field of view The field of view of our images was selected 
thanks to the following relation: 


FOV = d max /£ min (A.3) 

Here, the smallest available baseline was B m i n = 32 m and 
dmax = 665 nm, hence a field of view of 4.2 mas. 


7 Available at http://self-cal.oca.eu 
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